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ON THE EXISTENCE OF FLAT ORTHOGONAL 

MATRICES 

PHILIPPE JAMING AND MATE MATOLCSI 


Abstract. In this note we investigate the existence of flat or¬ 
thogonal matrices, i.e. real orthogonal matrices with all entries 
having absolute value close to ^=. Entries of ±-^= correspond to 
Hadamard matrices, so the question of existence of flat orthogonal 
matrices can be viewed as a relaxation of the Hadamard problem. 


1. Introduction 

Let M be a real orthogonal matrix of size n x n. We are interested 
in the smallest and largest modulus among the entries of M: 

(1) Im min and um max 

1 <i,j<n 1 <i,j<n 

and, more precisely, in estimating the maximal possible value for Im 
and the minimal possible value of %. In other words, we want to 
estimate the following quantities: 

(2) l n := max Im and u n := min um , 

MeO(n) MeO(n) 

where 0[n) denotes the group of orthogonal matrices. 

In order to determine l n (respectively, u n ) one has to control the 
lowest (resp. uppermost) absolute value in an orthogonal matrix M. 
It is also natural to ask whether one can control both quantities simul¬ 
taneously. For this purpose we introduce the measure of “flatness” of 
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an n x n orthogonal matrix M as 

f M = minis > 0 : —=- < I,j < ~t_ £ , for all 1 < i, j < n}, 

\/n V'n 

and we are interested in how flat an orthogonal matrix can be: 

(4) f n ■= min f M . 

MeO(n) 

ft is trivial to see that L < -4=, u„ > -4= for all n, and the famous 
Hadamard conjecture states that if 4|n then l n = u n = It is 

also easy to see that f n > unless there exists an n x n Hadamard 
matrix (Proposition 12. ll below). This note was motivated by the natural 
question of L. Baratchart (personal communication), as to whether 
l n > f° r some absolute constant c. The answer to this question 
is positive and provided by Theorem 12.21 below. We will also prove 
n„ < ^ with another absolute constant c in Theorem 12.31 However, 
we list here some more restrictive questions on flat orthogonal matrices, 
all of which we have only partial answers to. 

Problem 1.1. Is it true that l n = (1 + o( 1))-^? 

Problem 1.2. Is it true that u n — (1 + o( 1))^? 

Can we control l n and u n simultaneously? 

Problem 1.3. Is it true that f n = o(l) ? 

These questions can be seen as relaxations of the Hadamard prob¬ 
lem, and the matrices leading to such bounds could be called almost 
Hadamard matrices. However, a different notion of almost Hadamard 
matrices was already introduced and considered in [3j [2] (in those pa¬ 
pers the emphasis is on various matrix norms and not on the entries 
of the matrix). Therefore, to make a clear distinction, we prefer to use 
the terminology flat orthogonal matrices here. 

Finally, circulant matrices play a special role in applications, there¬ 
fore we can add this as an additional constraint. 

Problem 1.4. Let ff irc := mm Me0 (n)nCirc{n) I'm, where Circ(n ) de¬ 
notes the set ofnxn circulant matrices and f M was defined m (J3J) . Is 
it true that ff irc = o( 1) ? 

We firmly believe that the answer to the first three problems is posi¬ 
tive, while we are undecided as to the last one. Recall that a conjecture 
of R.yser asserts that there are no n x n circulant Hadamard matrices 
if 77 , > 4. 
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The remaining of this paper is split into two sections. The first one is 
devoted to general constructions that lead to bounds on u n and l n valid 
in arbitrary dimension. We then devote the last section to improved 
bounds when the size of the matrix have various arithmetic properties. 


2. General constructions 


We begin by a simple result which shows that the flatness parameter 
f n cannot be expected to be very small. In other words, requiring f n 
to be very small is equivalent to requiring a Hadamard matrix of order 
n to exist. 


Proposition 2.1. Let n > 3 and let e > 0 be such that 

1 jn if n is odd 


£ < 


2/n if n is even 


Assume there exists an orthogonal matrix M such that for every j, k = 


( 5 ) 


1 — £ 


n 


1/2 


< \m jtk \ < 


l+£ 


n 


1/2 


Then n is a multiple of 4 and there exists a Hadamard matrix of order 

n. 


Proof. Let S denote the matrix defined by the sign of the entries of M, 
i.e. Sj t k = ^signm^fc. Consider two rows Sj and s fc of S. Let r denote 
the number of columns where the entries in Sj and s^. match, and n — r 
where they differ. Then, for the corresponding rows rrq, m*, of M we 
have 


0 


(m i? m k ) 

{<-(1 + ^) 

< n 

^ n 


n — r 


n 

n — r 


n 


(1-e) 
(1 + e) 


2 r — n + en 


2 r — n — en 
n 


However, if n is odd the interval (2r — n — en, 2 r — n + en) does not 
contain zero (or any even number), a contradiction. If n is even, the 
interval (2r — n — en, 2 r — n + en) contains zero if and only if r — 
in which case the corresponding rows s j, s k are also orthogonal, and we 
conclude that S' is a Hadamard matrix. □ 


We continue with a simple block construction which proves that l n 
is at least as large as 

Theorem 2.2. For any dimension n we have l n > -Aj=. 
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Proof. Assume n = 2 r + q where q < 2 r , and introduce the notation 
s = 2 r — q for brevity. Let H be a Hadamard matrix of order 2 r (of 
course, such a matrix exists, e.g. the tensorial power F 2 0r ). We use 
the normalization that H has entries ±-4= instead of ±1. Let M be 

V 2 r 

the extension of H by an identity matrix of order q in the lower right. 
Split M into blocks of size s and q and q as follows (the indexes simply 
indicating the sizes of the blocks): 



" H S)S 

H s . q 

0 

M — 

Hq, s 

H q , q 

0 


0 

0 

h 


Let U denote the following orthogonal block-matrix: 



r is 

0 

0 ■ 

u = 

0 

L / 

LJ 


0 

V 2 1 

_L / 

V2 1 J 


and let M := U T MU. A direct calculation shows that 



Hs,S 

1 N 

V2 Us ’1 

1 N 

V2 Us ’ q 

M — 

1 TT 

2<A» + 4) 

+ Q 


I rr 

3 + h) 

3 (H,., + /,) 


The smallest modulus among the entries of M is which 

proves the theorem. (Note, however, that the largest appearing modu¬ 
lus is \ + y=, so this construction gives no indication with respect to 
Problem 11.31 1 □ 

Proving that u n < is also fairly easy, as one can use a block- 
diagonal construction. 

Theorem 2.3. For any dimension n we have u n < . 

Proof. If A 1 G 0(ni) and A 2 G 0(n 2 ) are orthogonal matrices then 
the block-diagonal matrix A\ © A 2 is an orthogonal matrix of order 
ni + n 2 . This implies u n+m < ma x{u n ,u m }. Also, Propositions 13.11 
and 13.21 below show that u p = l4 ^ 1 ' 1 whenever the dimension p is a 
prime. For general n we must invoke the following weak-Goldbach 
type result from number theory: for every e > 0, every large enough 
odd number n can be written as a sum of three primes n = pi +p 2 +p% 
where each pi lies in the interval [(1 — e)nj 3, (1 + e)nj 3], while every 
large enough even number n can be written as a sum of four primes n = 
P 1 +P 2 +P 3 +P 4 where eachp,; lies in the interval [(1— £)n/4, {l+efn/A\. 
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This follows from the method of Vinogradov |8j. Hence, the block- 
diagonal construction together with this weak-Goldbach type result 
implies u n < if n is odd, and u n < 2H ^i 1) if n is even. □ 

Remark 2.4. In connection with l n and u n it is natural to examine 
how the entries of a typical random orthogonal matrix behave. It is 
well known that for any coordinate z of a random unit vector we have 
Pr(\z\ > < e _< 7/2 (this can be seen by upper bounding the area of a 

spherical cap of radius r by that of a sphere of radius r). This implies, 
via a simple calculation, that um < c v^g n f or a random orthogonal 
matrix M, with high probability. Therefore, a random orthogonal ma¬ 
trix is typically “not far” from the bound ^ given in Theorem 

12.31 On the contrary, the lowest absolute value in a random unit vector 
is smaller than - with high probability (this can be seen by generat¬ 
ing a random unit vector as the normalized vector of n independent 
Gaussian random variables). The same holds, a fortiori, for a random 
orthogonal matrix, which shows that the lowest entry is typically very 
far from the optimal bound given in Theorem 12.21 □ 

Remark 2.5. If we assume that the Hadamard conjecture holds then 
In = u n = -^= for all n divisible by 4. By the simple construction of 
Proposition 13.1 1 we also conclude that for n = 3 ( mod 4) the quantities 
l n and u n are both of the magnitude (1 + o(l))^. However, we could 
not prove such a statement for the case n = 1, 2 ( mod 4). 

3. Specific constructions 

In the rest of this note we give some positive partial results with 
respect to Problems 11.1111.3111.41 As noted in (3J Section 3] symmetric 
balanced incomplete block designs give rise to orthogonal matrices with 
two entries. If the parameters of the block design are suitable then the 
entries will be close to ±^=. Namely, we have the following special 
case. 

Proposition 3.1. (i) If the dimension n is such that a Hadamard 
matrix H of size (n+1) x (n + 1) exists, then there exists an orthogonal 
matrix M of size n x n with all entries having modulus (1 + o(l))^=. 
In particular, this is the case if n = p r where p = 4k — 1 is a prime 
and r is odd. 

(ii) If the dimension n is a prime of the form 4k —1, then M can be 
chosen to be circulant. 

Proof, (i) We can assume that the Hadamard matrix H is in standard 
form, i.e. the first row and column of H consist of l’s. Delete the 
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first row and column of H , and in the remaining matrix H' replace the 
entries +1 and —1 by the variables x and y, respectively. Orthogonality 
of any two rows of H' is now equivalent to r ^x 2 + r ^-y 2 + xy = 0, 
and the unit length of the rows of H' is ensured by x 2 + ,J -^-y 2 = 1. 
This system of equations admits the (non-unique) solution 

-1 

x = . - 

\Jn + 1 — 1 

2 1 
^ \Jn + 1 \Jn — 1 — 1 

and we define M as the n x n matrix with these values. Note that 
both x and y have the order of magnitude (1 + o(l))-^. In fact, the 
error term o(l) here has the order of magnitude We remark 

that M corresponds to the Hadamard design associated to H , and this 
construction is a special case of the one described in [3, Section 3]. If 
n = p r where p = 4k — 1 is a prime and r is odd, then a Hadamard 
matrix of size (n + 1) x (n + 1) exists by the Paley construction. 

(ii) If n is a prime of the form 4k — 1 then the Paley construction 
essentially leads to a circulant matrix. Namely, consider the following 
circulant matrix M corresponding to the quadratic character of F p : 
[M] t j = x if i — j is a quadratic residue, and [M]^ = y if i — j is a 
non-residue or zero. The above values of x, y ensure orthogonality of 
M. □ 

When n is a prime of the form 4k + 1, the construction is much less 
trivial, as described below. 

Proposition 3.2. If the dimension is a prime p — 4k + 1 7 then there 
exists a circulant orthogonal matrix M of size p x p with all entries 
having modulus (1 + o(l))-^=. 

Proof. For this proof it will be convenient to identify the cyclic group 
with the numbers ( ~ p 0 +1 ,..., ^p-), and use coordinates of vectors 
accordingly. Also, let Q,NQ C 7L V denote the set of quadratic residues 
and non-residues, respectively (0 is not included in either Q or NQ). 
Consider the vector v = (v - P +i ,..., v P +i ) given by the quadratic char¬ 
acter, i.e. Vo = 0, and for nonzero j we have Vj = ±1 according to 
whether j G Q or j G NQ. Note that v is symmetric, Vj = v_j, be¬ 
cause p — 4k + 1. Note also that v(fc) = Vje 2mjk ^ p = ±^/p for 
k 0, and v(0) = 0. We will need a random modification of v. The 
construction is analogous to the one given in [6j Theorem 9.2], 
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Fix e > 0, and let p = A=. As NQ is symmetric, it can be written as 
a disjoint union of a set H and its negative, NQ — HU —H, in other 
words NQ = U ye n{y, ~y}- For each y e H consider independent 
random variables £ y such that P(£ y = 1/2) — p, P(£ J/ = 0) = 1 — p. 
For y = 0 set P(£ y = 1/4) = p, P(£ y = 0) = 1 — p. Consider the 
random vector w = (w -p+i , ..., w P -i ) given by w y = ^ y if y 6 H or 
y = 0, w y = C,- y if —y e H , and w y — 0 if y e Q. Let us evaluate the 
Fourier transform of the random vector w. 


(6) E(w(fc)) 

and, for all k, 

(7) D 2 (wl») = 




E 1 

yeNQ 


0 2jriyk/p 


y/p /4 if k — 0, 
±1/4 if k ± 0 



g 27r iyk/p _|_ e ~2niyk/p\2 j < - 


A <~ 




Vp- 


We invoke here a large deviation inequality of Chernov as stated 
in (?) Theorem 1.8]: let A/,..., X n be independent random variables 
satisfying \X t — E(Xj)| < 1 for all i. Put A" = X 1 + ... + X n and let 
a 2 be the variance of X. For any t > 0 we have 

P(|X - E(AT)| >ta)< 2 max (V* 2/4 , e~ tal2 ^ . 


Using this estimate with t = p £ we obtain 


(8) P (jw(0) - y/p/A\ > \p li+£ ^j < 2e ?4 , 

(9) P ^|w(fc) || > < 2e _E ^ for all k ^ 0. 

Therefore, with high probability none of the above events occur, and 
we have |w(0) — yjp/ 4 | < 0(p* +£ ) and |w(fc)| < 0(p 4 +£ ). Fix such 
a favourable vector w. We can also assume without loss of generality 
that £ 0 = 1/4 and hence w 0 = 1/4 (we are free to change w 0 from 0 to 
1/4, if necessary, without altering the order of magnitude of w). 

Finally, consider the vector z = v ± 4w, and let u = lz. The vector 
z is unimodular (some of the —1 entries in v were changed to ±1 and 
the value at 0 was changed to ±1). Therefore the circulant matrix M 
with first row u is orthogonal. Also, the entries of u are all of absolute 
value -F(l ±p - 4+ e ) ; by construction. 
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We remark that a similar construction works if the dimension p is of 
the form 4k — 1 but the result is inferior to Proposition 13 .1 1 in the sense 
that the error term is larger. □ 

By combining the results of the propositions above we can answer 
Problem 11.31 and 11.41 for dimensions n which are composed of large 
prime factors. 

Corollary 3.3. For any fixed m let N m denote the set of positive in¬ 
tegers n = 2 s pi.. ,p r such that each odd prime factor pj > n 1//m (the 
primes may appear with multiplicity). 

(i) For every m and every n G N m there exists an orthogonal matrix 
M of size n x n with all entries having modulus (1 + o m (l))^. 

(ii) If all the odd primes pj appearing in the factorization of n are dis¬ 
tinct and s = 0 or s = 2 then there exists a circulant orthogonal matrix 
M of size n x n with all entries having modulus (1 + o m (l))^. 

Proof, (i) Let n = 2 s pi.. ,p r e N m . Then r < m by the definition of 
N m . Let H denote a Hadamard matrix of size 2 s (such matrix exists, of 
course, e.g. H = F ® s ), and let M Pj denote the matrices corresponding 
to the primes pj, as constructed in Propositions 13.II and 13.21 Let M be 
the tensorial product of all the matrices H and M p . (where M p . is taken 
with the same multiplicity as pj in n ). As r <m for any n e N m , the 
errors do not accumulate, and we still have |[M]jj| = (1 + o m (l))^. 

(ii) The matrices M p . are circulant by the constructions of Propo¬ 
sitions 13.11 and 13.21 Let Xj denote the hrst row of M Pj . If s = 0, 
n = pi ... p r with all pj distinct, then 7L n = Z pi ...Z Pr . Let k = 
(fci,..., Ay) e Z pi ... Z Pr and let y(k) = n^i x j(%)- Then y gener¬ 
ates a circulant matrix which is orthogonal (because all x.,■ are uni- 
modular, and hence so is y), and the entries of y are of absolute value 
(1 + o m (l))^=. When s = 2 we can incorporate the 4x4 real circulant 
Hadamard matrix in the same manner. □ 

Remark 3.4. Problem 11.41 concerning the circulant case has an in¬ 
teresting connection to ultraflat polynomials. It is well-known that 
x = (x \,..., x n ) generates a circulant orthogonal matrix if and only 
if the Fourier transform x = (wi,... ,w n ) is unimodular on Z n , i.e. 

| vjj | = |iDfc| for all j,k. If x is real, then w is conjugate symmet¬ 
ric, i.e. Wj = w n _j. We also want that all \xf ~ Considering 
wi,...,w n as variables we are led to the problem of constructing a 
polynomial P(z) = WjZ 3 where Wj = w n _j and | w 3 \ = 1 (after 

re-normalization), such that P(z) is “flat” at the nth roots of unity, 
i.e. \P(uj j )\ = (1 + o(l))i/n, where oj = e 2m i n and j = 0,.. .n — 1. 
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Dropping the restriction Wj = w n -j one can require P(z) to be flat 
all over the unit circle (i.e. \P(z)\ = (1 + o(l))y/n for all \z\ = 1), 
and such polynomials are called ultraflat. The existence of ultraflat 
polynomials was proven by Kahane [5jj. However, the extra condition 
Wj = w n -j prevents P(z) from being ultraflat as shown by Remark 5.1 
in [4]: the restriction Wj = w n -j implies maxy-ei \P(z)\ > (1 + e)y/n 
with e = -y/4/3 — 1. However, this does not answer Problem 11.41 be- 
cause we require P(z) to be flat only at the nth roots of unity. Problem 
P can therefore be regarded as a discretized version (at the nth roots 
of unity) of the question of existence of ultraflat polynomials with the 
self-conjugacy restriction Wj = w n -j- 

Notice also that Corollary 13.31 gives an affirmative answer to Prob¬ 
lem 11.41 for some dimensions n (e.g. when n = 4 p, p being a prime). 
Therefore, trying to prove Ryser’s conjecture on the non-existence of 
circulant Hadamard matrices by giving a negative answer to Problem 
P cannot possibly work. In the other direction, hoping to construct 
a circulant Hadamard matrix by first constructing a circulant flat or¬ 
thogonal matrix and then modifying its entries is also rather naive and 
hopeless in our opinion. □ 

We saw in Proposition 13.11 that given any Hadamard matrix H one 
can reduce the dimension by 1, and construct a flat orthogonal matrix. 
It is natural to try also to increase the dimension by 1. The general 
construction given in Proposition 12. H allows to do so, but will introduce 
an entry of size 1/2, thus destroying the flatness. Our last specific 
construction concerns the increase of dimension by 1 without destroying 
flatness, but it only works under some restrition on the dimension. We 
recall that a Hadamard matrix is called regular if the row sums and 
column sums of H are all equal. 

Proposition 3.5. Assume the dimension n is such that a regular 
Hadamard matrix H of order n exists (this implies n = 4 k 2 for some 
k). Then there exists an orthogonal matrix M of size (n + 1) x (n + 1) 
with all entries having modulus (1 + o(l)) J~ +1 • 

Proof. Let b = —— and replace the positive entries of H by — 

b, while the negative entries with —A= — b (the sign of b is purposefully 
negative in both cases). Next, extend this matrix by a new row and 
column filled with entries a = ^=, and let M be the arising matrix. 
The entries of M are of modulus (1 + o(l))^== (in fact, the error 
term is of the order ^), and an easy calculation shows that M is 
orthogonal. □ 
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We end this note by emphasizing that Problems 11.11 11.31 11.41 all 
remain open for dimensions n with small prime factors, i.e. for dimen¬ 
sions n which are not covered by Corollary 13.31 
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